Nonlocal response in plasmonic waveguiding with extreme light confinement 
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We present a novel wave equation for linearized plasmonic response, obtained by combining the 
coupled real-space differential equations for the electric field and current density. Nonlocal dynam- 
ics are fully accounted for, and the formulation is very well suited for numerical implementation, 
allowing us to study waveguides with subnanometer cross-sections exhibiting extreme light confine- 
ment. We show that groove and wedge waveguides have a fundamental lower limit in their mode 
confinement, only captured by the nonlocal theory. The limitation translates into an upper limit for 
the corresponding Purcell factors, and thus has important implications for quantum plasmonics. 

PACS numbers: 78.67.Uh, 78.67.Lt, 71.45.Lr, 73.20.Mf, 41.20.Jb 
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Wave propagation along dielectric waveguide struc- 
tures has over the years been extended also to plas- 
monic systems with waveguide modes in the form of 
surface-plasmon polaritons. Plasmonic waveguides have 
attracted considerable attention during the past decade, 
primarily due to their ability to support extremely con- 
fined modes, i.e., modes that do not exhibit a diffraction- 
limited cutoff for progressively smaller waveguide cross 
sections but transform themselves into their electrostatic 
counterparts [1]. Investigations of nanowire [2], groove [3] 
and wedge [4] waveguides, shown to ensure extreme light 
confinement, raise a natural interest in the influence of 
nonlocal effects on strongly confined plasmonic modes [HI . 
Waveguiding by metal nanowires tej and more recently 
plasmonic focusing by conical tips [3, [H have been stud- 
ied in the context of nonlocal response. However, with 
the exception of few analytical studies of simple planar 
geometries [ol, [lo[ , nonlocal effects in the dispersion prop- 
erties of complex waveguides remain unexplored, a cir- 
cumstance that can partly be explained by the added 
complexity due to nonlocal effects as compared to the 
widespread framework of the local-response approxima- 
tion (LRA) [nj. 

There is also another good reason to look for non- 
local effects in extreme light confinement. Subwave- 
length mode confinement implies large effective Purcell 
factors and thereby strong coupling of single emitters 
to nearby plasmonic waveguide modes [l2|. The latter 



opens a doorway to quantum optics with surface plas- 
mons, including the possibilities for realization of single- 
photon transistors [l3| and long-distance entanglement of 
qubits [l^. Since one would expect that the plasmonic 
mode confinement is fundamentally limited by nonlocal 
effects, similarly to nonlocal limits in the field enhance- 
ment of localized plasmon excitations [HI, [l6| , studies of 
the plasmonic mode confinement beyond the LRA are of 



great interest for quantum plasmonics. 

In this Letter, we derive a novel wave equation which 
fully takes into account the nonlocal dynamics of an 
often-employed hydrodynamical model (HDM). We ap- 
ply the wave equation to plasmonic waveguides (Fig. [1]) 
with extreme light confinement, defined by the sub- 
nanometer dimensions of the waveguide cross section. 
After stringent bench-marking of our approach against 
the analytically tractable case of nanowires with circu- 
lar cross-section, we analyze in detail groove and wedge 
waveguides and demonstrate the existence of fundamen- 
tal limits in their mode confinement and Purcell factors, 
imposed by the nonlocal effects. 

The nonlocal response, or spatial dispersion, is a con- 
sequence of the quantum many-body properties of the 
electron gas, which we here take into account within a 
semi-classical model 17-2o|. In this model the equation- 
of-motion for an electron in an electrical field is supple- 
mented with a hydrodynamic pressure term originating 




FIG. 1: Generic plasmonic waveguiding geometries with wave 
propagation in the 2;-direction and extreme transverse con- 
finement in the xy-plane due to subnanometer geometric di- 
mensions, e.g. the nanowire radius R or the edge radius-of- 
curvature r. 
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from the quantum kinetics of the electron gas. By hn- 
earization, the plasmonic response is governed by the 
fohowing pair of coupled real-space differential equa- 



tions 



21|: 



Vx VxE(r)== (f)'E(r) + iwAioJ(r), (la) 

;;;(^V[V-J(r)]+J(r)=<T(r)E(r). (lb) 

Here, the term V[V-J] = VxVxJ + J is a correc- 
tion to Ohm's law and scales as = {3/5)vp within the 
Thomas-Fermi model [22|] with vf being the Fermi veloc- 
ity. For simplicity we neglect here any interband effects 
present in real metals; these can be included straightfor- 



wardly [23|, • numerical solutions we will con- 

sider Drude parameters appropriate for silver [25|. As- 
suming a hard-wall confinement associated with a high 
work function, the boundary conditions for the current 
at the metal surface become particularly simple: the tan- 
gential component is unrestricted while the normal com- 
ponent vanishes due to the current continuity and vanish- 
ing of all electron wave functions at the surface lol, l21 |. 

For analytical progress one can eliminate the current 
from Eq. ([Taj), thereby arriving at an integral equation 
where a dyadic Green's function accounts for the nonlocal 
dynamics of the electron gas [26| . Alternatively, the cou- 
pled equations ([Taj) and ([Tb|) form a natural starting point 
for a numerical treatment of arbitrarily shaped metal- 
lic nanostructures, e.g., with a state-of-the-art finite- 
element method 
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27j . Recently, we employed this ap- 
proach to study field enhancement and SERS in groove 
structures [l5|. However, for waveguiding geometries we 
seek solutions of the form E(r) ex ex.p{ikzz) leading to 
an eigenvalue problem for kz{uj) with a six-component 
eigenvector {E, J}. In that context the coupled-equation 
formulation is numerically less attractive. Here, instead, 
we eliminate the current from Eq. (|lb[) . a procedure that, 
after straightforward manipulations using standard vec- 
tor calculus 0, results in an appealingly compact, but 
yet entirely general nonlocal wave equation: 

VxVxE(r) = (^)'e,,(r)E(r), 



(2a) 



.(r)=e„(r) + 
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uj{uj-\-i/ t) 



(2b) 



Here, the operator £NL(r) contains the nonlocal effects. In 
the limit /3 ^ 0, £NL(r) reduces to the usual Drude dielec- 
tric function ^^(r) = l^i(j{Y) / {eqUo) = 1 — UJp{T)/[ljj{uJ -\- 

i/r)] used in the LRA. Thus, with a simple rewriting we 
have turned the coupled- wave equations into a form rem- 
iniscent of the usual wave equation, with all aspects of 
nonlocal response contained in the Laplacian term /3^V^ 
in £nl(i*). This is the main theoretical result of this Let- 
ter. Clearly, the single-line form is beneficial for the con- 
ceptual understanding and further analytical progress. 
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FIG. 2: Fundamental waveguide mode of a cylindrical silver 
nanowire embedded in air. (a) Dispersion relation uj{kz) and 
(b) normalized effective mode area within the HDM for the 
nanowire radius R = 2nm (blue) and 4nm (green), respec- 
tively, showing excellent agreement between numerical solu- 
tions of Eq. (|2]) (solid points) and analytical results (solid 
lines). For comparison, the red-dashed curve shows the uni- 
versal result of the nonretarded LRA, with its large- /c^ limit- 
ing value of ujp/\/2 indicated in (a) by the horizontal line, (c) 
Radial distribution of the electric field |Ep | at cj = 0.6 cjp for 
R = 4nm, contrasting the continuous field variation in the 
HDM with its usual boundary discontinuity in the LRA. 



as well as for numerical implementations: the additional 
Laplacian does not add any complications beyond those 
already posted by the double-curl operator on the left- 
hand side equation. The result in Eq. ([2j) is instrumental 
for efficient and accurate numerical eigenvalue solutions 
in waveguiding geometries with arbitrarily shaped waveg- 
uide cross sections. 

We now apply the developed formalism to the waveg- 
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uide configurations of Fig. [T] wiiicii can provide extreme 
ligtit confinement i) metal nanowires witii circular 
cross sections [H where analytical solutions 0] are avail- 
able for benchmarking of the numerics, ii) grooves in 
metal |3i], and in) metal wedges [^]. In addition to the 
usual mode characteristics, effective index and propaga- 
tion length, we also evaluate the effective mode area: 
^eff = Veii/L^ where Kff is the effective mode volume 
associated with the Purcell effect, i.e., 



^eff 



dxdyu{T) 



maxy^ {^/(r)} 



(3) 



where u{t) is the electromagnetic energy functional [24 1. 
The cross-sectional integral extends over the volumes Vm 
and Va occupied by metal and air, respectively, while the 
evaluation of the maximal energy density is restricted to 
the air region where dipole emitters can be placed. 

The dispersion curves and effective mode areas (nor- 
malized to the nanowire cross section) calculated for sil- 
ver nanowires of different radii [Fig. [2ja,b)] exhibit a 
blueshift and increased mode area (for fixed kz) when 
taking nonlocal effects into account. They show excellent 
agreement with the corresponding analytical results 0. 
Importantly, nonlocal dynamics influences strongly the 
mode field distribution [see Fig. [2fc)], because, contrary 
to the LRA case, the normal component of the electrical 
field within the HDM is continuous across the interfaces 
(this is a special case for a Drude metal without inter- 
band effects and surrounded by vacuum [2^ ) . It is indeed 
seen [Fig. [2fc)] that |Ep| is discontinuous on the bound- 
ary in the local case, while it varies continuously across 
the boundary in the nonlocal case. This variation occurs 
in a region extending over ^0.1 nm, that is of the order 
of the Fermi wavelength of silver. 

The results for cylindrical nanowires, while demon- 
strating the main effects of nonlocal dynamics on the 
mode characteristics, indicate that the quantitative 
changes are modest even for very small radii (Fig. [2]). 
In order to explore fundamental limitations, one has to 
consider the limit of vanishing radii of curvature. While 
subnanometer radii appear unrealistic for nanowires, 
fabrication of grooves cut in metal and metal wedges, 
e.g., by nanoimprint lithography |28i], can in fact result 
in nm-sharp edges with corresponding nm-sized wedge 
modes [1]. We expect that nonlocal effects then come 
into play. 

Rather surprisingly, the mode effective index and prop- 
agation length calculated for silver grooves and wedges 
(Fig. [3]) exhibit even weaker influence of the nonlocal 
effects as compared to the case of nanowires (Fig. [2]). 
In fact, there is no noticeable difference between the 
LRA- and HDM-based results obtained for 1-nm-radius 
of edges. In the limit of mathematically sharp edges, the 
mode effective index becomes only slightly larger and the 
propagation length slightly smaller than those calculated 



for 1 nm edge radius (Fig. [3]). We explain this result by 
the fact that groove and wedge plasmonic modes are only 
partially affected by the very tip, being distributed also 
and predominantly over flat edges (see insets in Fig. [3]). 

The situation changes drastically when one considers 
the mode confinement, using the mode area associated 
with the Purcell factor, Eq. (j3j). We recall that the field 
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FIG. 3: Effective index (left axis) and propagation length 
(right axis) versus wavelength for the fundamental mode in 
complimentary (a) V-groove and (b) wedge silver waveguides, 
both with an opening angle of 30° . The nonlocal results (solid 
circular symbols) obtained with Eq. ([2]) are contrasted to the 
LRA (open circles), with dashed lines serving as eye guides. 
Results for mathematically sharp structures with r = (blue 
solid circles) are contrasted to finite rounding with r = 1 nm 
(red open circles). Insets show field- intensity distributions 
(white scale bars are 1 nm long) calculated within the HDM 
(A = 600 nm) for infinitely sharp edges. The fingerprint of 
nonlocal effects is clearly visible as the field penetrates into 
the metal by a distance of the order of the Fermi wavelength 
of silver. 
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FIG. 4: Normalized mode area versus wavelength for the fun- 
damental mode in complimentary (a) V-groove and (b) wedge 
silver waveguides, both with opening angles of 30°. The HDM 
results (sohd symbols) are contrasted to the LRA (open cir- 
cles) for r = 1 nm (red) and r = 0.2 nm (green). Results for 
mathematically sharp structures with r = (blue solid cir- 
cles) define a lower limit in the HDM (grey-shaded regions 
are inaccessible). For the LRA, the r = 0.1 nm results (ma- 
genta) exceed this limit and the mode area tends to zero when 
r ^ 0. Insets show field- intensity distributions (white scale 
bars are 5 nm long) at A = 600 nm. The LRA intensities are 
with rounding r = 1 nm, while r = is used for the HDM 
maps. 



enhancement calculated within the LRA grows without 
bound for progressively sharper pointed structures while 
it remains finite when calculated within HDM 
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Analogously, in the present case, one may expect that the 
mode area calculated within the LRA decreases without 
bound for a decreasing edge radius, while it may satu- 
rate within the HDM. LRA-based simulations for sub- 
nanometer radii of edges show (Fig. S]) that the mode 
area indeed tends to zero, without any apparent satura- 



tion. This trend is more pronounced for wedges than for 
grooves, because the wedge geometry ensures generally a 
better mode confinement [cf. Figs. HJ^a) and (b)], as was 
also noted previously At the same time, the simu- 
lations conducted within the HDM demonstrate clearly 
the existence of a lower bound for the mode area which 
remains finite even for mathematically sharp edges (blue 
circles in Fig. 2]). The associated Purcell factors can be 
estimated by inverse of the normalized mode areas dis- 
played in Fig. |4] [29^. Thus, our calculations show that 
there is a fundamental limit for the maximum Purcell 
factors achievable with plasmonic waveguides. It is inter- 
esting that the upper limit of Purcell factors evaluated 
in this way decreases noticeably in the long-wavelength 
regime. This feature is related to a general weakening of 
all plasmonic effects, including waveguiding [l], for longer 
wavelengths (with metals approaching the limiting case 
of perfect conductors). At the same time, in the case 
of wedges, these factors remain substantial even at tele- 
com wavelengths, with the propagation lengths becom- 
ing considerably long (Fig. [3]) and amenable for circuitry 
application. It should also be borne in mind that the 
plasmonic field confinement in both grooves and wedges 
increases for smaller opening angles [l, 0] ^ so that even 
larger Purcell factors can be achieved, albeit at the ex- 
pense of shorter propagation. 

In conclusion, using a novel wave equation accounting 
for nonlocal dynamics, we considered plasmonic waveg- 
uides with extreme light confinement and demonstrated 
the existence of a fundamental limit in their mode con- 
finement imposed by nonlocal effects. Our results imply 
fundamental limitations in the corresponding Purcell fac- 
tors, showing at the same time the possibility of achieving 
very high Purcell factors with V-groove and wedge waveg- 
uides that ensure sufficiently long propagation lengths for 
applications in quantum plasmonics. 
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